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Simulation  of  the  northern  Adriatic  circulation  Approved  for  Public  Release 
during  winter  2003  Distribution  Unlimited 
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[i|  Numerical  simulations  of  the  Adriatic  Sea  were  conducted  with  the  Navy  Coastal 
Ocean  Model  during  the  Adriatic  Circulation  Experiment  in  the  fall  and  winter  of 
2002/2003.  and  results  were  compared  with  observations.  The  ocean  model  used  a  1-km 
resolution  grid  over  the  entire  Adriatic  Sea.  Model  forcing  included  atmospheric  fluxes 
from  the  Coupled  Ocean-Atmosphere  Mesoscale  Prediction  System  (COAMPS),  tides, 
boundary  conditions  from  a  global  model,  and  freshwater  river  and  runoff  inflows.  Model 
tidal  elevation  showed  good  agreement  with  International  Hydrographic  Organization 
station  data,  and  model  tidal  currents  showed  good  agreement  with  tidal  currents 
determined  from  acoustic  Doppler  current  profiler  (ADCP)  measurements.  Detided  model 
currents  showed  good  agreement  with  ADCP  currents  with  rms  errors  along  the  principal 
variance  axes  ranging  from  6  to  12.9  cm/s  and  correlations  ranging  from  0.16  to  0.81. 

Correlations  between  model  and  ADCP  currents  along  the  minor  variance  axes  were 
generally  low.  Comparison  of  the  model-simulated  temperature  and  salinity  profiles 
during  January  and  February  with  conductivity-temperature-depth  measurements 
indicated  that  the  model  captured  some  of  the  spatial  structure  of  the  observed  fields.  The 
model  response  to  several  bora  wind  events  in  January  and  February  showed  a  recurring 
pattern  of  cyclonic  and  anticyclonic  gyres  that  generally  agreed  with  observations  and 
reflected  the  pattern  of  wind  stress  curl  from  the  COAMPS  wind  stress  forcing.  During 
strong  bora  forcing,  two  large  cyclonic  circulation  gyres  form  in  the  northern  Adriatic  with 
a  smaller  anticyclonic  circulation  between  them  near  the  Istrian  Peninsula.  Additionally,  a 
couple  of  large  meanders  frequently  occur  within  the  Eastern  Adriatic  Current  southeast  of 
Kvamer  Bay,  and  these  meanders  sometimes  close  to  form  small  gyres. 

Citation:  Martin.  P.  J..  J.  W.  Book,  and  J.  D.  Doyle  (2006).  Simulation  of  the  northern  Adriatic  circulation  during  winter  2003. 
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1.  Introduction 

[2]  There  has  been  a  long-standing  interest  in  theoretical 
and  modeling  studies  of  the  Adriatic  Sea  as  it  is  charac¬ 
terized  by  strong  wind  and  river  forcing  and  is  a  site  of 
dense  water  formation  for  the  eastern  Mediterranean  Sea 
[Artegiani  et  at,  1997a,  1997b;  Cushman-Roisin  et  at ,  20011. 

[y]  The  general  circulation  of  the  Adriatic  is  cyclonic 
with  southeastward  flows  along  the  western  side  of  the  sea 
and  northwestward  flows  along  the  eastern  side  [Orlic  et  at, 
1992],  Three  smaller  cyclonic  patterns  are  frequently  ob¬ 
served  within  the  overall  cyclonic  circulation  [Artegiani  et 
at ,  1997b;  Paula  in^  2001] 

[4]  In  winter,  the  Adriatic  is  subject  to  a  strong  wind 
called  a  bora,  which  is  a  cold  northeasterly  wind  that 

'Ckctfti  Dynamics  and  Prediction  Branch,  Oceanography  Division, 
Naval  Research  Laboratory,  Stearns  Space  Center,  Mississippi,  USA 

^Ocean  Sciences  Branch,  Oceanography  Division,  Naval  Research 
Laboratory,  S tennis  Space  Center,  Mississippi,  USA. 

'Atmospheric  Dynamics  and  Prediction  Branch,  Marine  Meteorology 
Division,  Naval  Research  Laboratory,  Monterey,  California,  USA. 

Tliis  paper  is  not  subject  to  U.S,  copyright. 

Published  in  2006  by  the  American  Geophysical  Union 


develops  in  the  lee  of  the  Dinaric  Alps.  A  number  of  studies 
suggest  that  the  bora  shares  some  common  characteristics 
with  downs  I  ope  windstorms  and  transcritical  hydraulic 
flows  [e.g..  Smith,  1 987;  Klemp  and  Durr  an,  1987;  Jiang 
ami  Doyle,  2005]  The  bora  winds  impact  both  the  water 
characteristics  and  the  circulation,  especially  in  the  shallow 
northern  end  of  the  Sea.  The  other  primary  wind  pattern  is 
the  sirocco,  which  blows  from  the  southeast  along  the  axis 
of  the  Adriatic. 

[5]  There  is  a  net  gain  of  fresh  water  within  the  Adriatic 
Sea,  unlike  the  Mediterranean  as  a  whole.  Total  mean 
freshwater  inflows  into  the  Adriatic  have  been  estimated  to 
be  about  5700  m  Vs  [Raicich,  1994].  The  largest  single  source 
of  fresh  water  in  the  Adriatic  is  the  Po  River,  which  has  an 
annual  mean  discharge  of  about  1500  m  Vs,  This  low-salinity 
water  flows  southeast  along  the  Italian  coast  [Hopkins  et  at, 
1999]  and  joins  a  current  turning  to  the  southeast  along  the 
western  Adriatic  slope  [Zore-Armanda  ei  at,  1996]  to  form 
the  Western  Adriatic  Current  (WAC).  A  less  intense  flow 
along  the  eastern  side  of  the  Adriatic  forms  the  Eastern 
Adriatic  Current  (EAC). 

fr]  Past  theoretical  studies  have  been  successful  at  de¬ 
scribing  the  bulk  features  of  the  oceanography  of  the 
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Adriatic  [e.g.,  Hendershott  and  Rizzoli,  1976;  Ortic  et  ai , 
1994].  However,  recent  numerical  modeling  studies  have 
begun  to  focus  on  the  more  complicated  details  of  the 
Adriatic,  recognizing  the  importance  of  the  fine  scales  of 
wind  forcing  [Pullen  et  ai ,  2003]  and  using  complex 
topography  and  models  [Beg  Paklar  et  ai ,  2001;  Zavatarelli 
et  ai,  2002]  to  simulate  more  realistic  conditions  and  study 
the  prevailing  dynamical  details.  Comparisons  with  data  for 
these  studies  have  been  limited  as  most  of  the  observations 
were  confined  to  single  point  samples  and  were  limited  in 
spatial  and  temporal  extent. 

[?]  An  extensive  international  collaborative  study  of  the 
northern  and  central  Adriatic  was  conducted  during  the  fail 
and  winter  of  2002/2003,  which  included  a  number  of 
observational  programs  and  numerical  modeling  efforts. 
Measurements  that  were  made  include  acoustic  Doppler 
current  profilers  (ADCPs),  conductivity -temperature-depth 
(CTD)  casts,  surface  drifter  measurements,  high-frequency 
radar  measurements,  meteorological  platform  and  buoy 
measurements,  towed  body  measurements,  and  other  cur¬ 
rent  mooring  measurements  [Lee  et  ai ,  2005].  Numerical 
modeling  efforts  included  high-resolution  meteorological, 
ocean,  wave,  and  sediment  model  simulations. 

[hJ  This  paper  describes  a  numerical  simulation  of  the 
Adriatic  that  was  conducted  with  the  Navy  Coastal  Ocean 
Model  tor  this  time  period  and  a  comparison  of  the  results 
from  the  simulation  with  some  of  the  ADCP  and  CTD 
observations.  The  main  focus  is  the  period  of  January- 
February  2003,  during  which  several  bora  events  occurred. 
The  following  sections  include  descriptions  of  (2)  the  ocean 
model,  (3)  the  ocean  model  setup,  (4)  the  observations  used 
for  validation  of  the  model  simulation,  (5)  the  model  results 
and  comparison  with  observations,  and  (6)  summary  and 
conclusions. 

2,  Ocean  Model 

M  The  ocean  model  used  here  is  the  Navy  Coastal  Ocean 
Model  (NCOM)  as  described  by  Martin  [2000],  with  some 
improvements  as  described  by  Morey  et  ai  [2003]  and 
Barron  et  ai  [2006].  This  model  is  similar  in  its  physics  and 
numerics  to  the  Princeton  Ocean  Model  (POM)  [Blum berg 
and  Mellor „  1987],  but  uses  an  implicit  treatment  of  the  free 
surface  and  a  hybrid  vertical  grid  with  sigma  coordinates  in 
the  upper  layers  and  (optionally)  level  coordinates  belowr  a 
user-specified  depth.  Only  the  sigma  coordinate  part  of  the 
grid  moves  with  the  free  surface, 

[10]  The  model  equations  include  a  source  term  that  can 
be  used  for  river  inflow's.  There  are  options  for  higher-order 
treatment  of  some  terms,  for  example,  third-order  upwind 
for  advection  [Holland  et  ai ,  1998],  which  was  used  for  the 
simulations  conducted  here.  Vertical  mixing  is  computed 
using  the  Mellor- Yanriada  Level  2  scheme  [Mellor  and 
Yarnada ,  1974],  which  is  modified  for  use  over  the  entire 
water  column  .  The  equation  of  state  is  that  of  Mellor  [  1991  ] 
as  used  in  POM. 

3.  Ocean  Model  Setup 
3X  Ocean  Model  Domain 

[11]  The  ocean  model  domain  consists  of  the  entire 
Adriatic  Sea  and  includes  the  Strait  of  Otranto  and  a  small 


part  of  the  northern  Ionian  Sea.  Three  different  computa¬ 
tional  grids  have  been  used  for  the  Adriatic  modeling:  a 
3-km  resolution  grid  with  a  1-km  nest  in  the  nonhem 
Adriatic  and  single  2 -km  and  l -km  resolution  grids  over 
the  entire  Adriatic.  Comparison  of  results  on  the  different 
grids  have  provided  some  indication  of  the  sensitivity  of  the 
results  to  the  model  grid  and  bathymetry  Results  presented 
here  are  from  the  1-km  grid  covering  the  full  Adriatic. 
Figure  1  show's  the  bathymetry  for  the  entire  domain 
covered  by  this  grid  and  Figure  2  shows  the  bathymetry 
in  the  northern  Adriatic  in  more  detail.  The  bathymetry  was 
derived  from  a  database  developed  by  the  Naval  Oceano¬ 
graphic  Office  and  nautical  chart  soundings  and  the  land-sea 
boundary'  was  derived  from  the  Generic  Mapping  Tools 
(GMT)  vector  shoreline. 

[12]  The  vertical  coordinate  used  with  the  1-km  grid 
consists  of  32  total  layers,  with  22  sigma  layers  used  from 
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Figure  2,  Model  bathymetry  in  nonhem  Adriatic  (contour 
interval  is  10  m)  with  mooring  locations  (labeled  open 
circles)  and  locations  oTCTDs  taken  in  January- February 
2003  (small  dots). 

the  surface  down  to  a  depth  of  291  m  and  level  coordinates 
used  below  291  m.  Hence  the  grid  is  like  a  regular  sigma 
coordinate  grid  in  water  shallower  than  291  m  and  is  similar 
to  a  level  grid  in  deeper  water.  The  vertical  grid  uses  a 
uniform  stretching  with  a  maximum  thickness  of  the  upper 
layer  of  2  m  and  a  maximum  depth  of  1262  m.  The 
minimum  bottom  depth  was  set  to  2  m, 

3*2.  Initial  and  Boundary  Conditions 

[n]  Initial  conditions  (1C)  for  sea  surface  height  (SSH), 
velocity,  temperature  (T),  and  salinity  (S)  were  interpolated 
from  a  hindcast  of  the  Naval  Research  Laboratory  (NRL) 
global  NCOM  model  \Barron  et  aLy  2004]  for  00Z, 
I  September  2002.  Boundary  conditions  (BC)  at  the  open 
boundary  in  the  northern  Ionian  Sea  were  provided  by  daily 
values  of  SSI  L  velocity,  T,  and  S  from  the  global  model. 
The  numerical  treatment  of  the  BC  includes  the  Flather 
radiative  BC  [Flather  and  Proctor ,  1983]  for  the  SSH  and 
depth-averaged  normal  velocity,  Orlanski  radiation  condi¬ 
tions  [  Or  Ians  ki,  1976]  for  the  tangential  velocities,  T,  and  S, 
and  a  relaxation  to  the  T  and  S  fields  of  the  global  model 
near  the  open  boundary.  The  normal  baroclinic  velocity  at 
the  open  boundary  is  computed  using  the  full  model 
equation  except  that  advection  is  only  computed  normal 
to  the  boundary  using  a  first-order  upwind  scheme. 

3.3.  Tidal  Forcing 

[14]  Tidal  forcing  was  provided  using  tidal  SSH  and 
depth-averaged  normal  and  tangential  velocities  at  the  open 
boundaries  from  the  Oregon  State  University  (OSU)  tidal 
databases,  which  are  derived  from  satellite  altimetry  data 
[Egbert  and  Erofeeva *  2003],  The  tidal  data  are  linearly 
added  to  the  BC  from  global  NCOM  (which  does  not 


include  tides).  Data  from  the  OSU  Mediterranean  tidal 
database  were  used  for  the  Kl,  01,  M2,  and  S2  constituents 
and  data  from  the  OSU  global  database  were  used  for  PI, 
QL  K2,  and  N2.  Tidal  potential  forcing  for  these  eight 
constituents  was  used  in  the  interior  of  the  model  domain. 

3.4.  Atmospheric  Forcing 

[is]  Atmospheric  forcing  was  obtained  from  the  atmo¬ 
spheric  component  of  the  Coupled  Ocean- Atmosphere  Me- 
soscale  Prediction  System  (COAMPS™)  [Hodur,  1997]. 
The  COAMPS  setup  for  the  Adriatic  consists  of  a  triply 
nested  grid  with  resolutions  of  36,  12,  and  4  km  [Pullen  et 
al ,  2003],  The  outer  grid  of  this  nested  grid  system  covers 
most  of  Europe  and  the  Mediterranean  and  the  inner  4-km 
grid  covers  the  entire  Adriatic  and  part  of  the  Tyrhennian 
Sea.  COAMPS  itself  is  nested  within  the  Navy  Opera¬ 
tional  Global  Atmospheric  Prediction  System  (NOGAPS) 
[Rosmond  et  al .,  2002]. 

[  16]  Atmospheric  forcing  was  provided  by  hourly  fields 
of  surface  air  pressure,  wind  stress,  solar  radiation,  net 
longwave  radiation,  and  precipitation  from  the  COAMPS 
4-km  grid.  Latent  and  sensible  heat  fluxes  were  computed 
with  standard  bulk  formulas  using  the  COAMPS  10-m  wind 
speed  and  2-m  air  temperature  and  humidity  and  the  ocean 
model  SST.  The  stability-dependent  Kondo  [1975]  drag 
coefficient  was  used  for  the  hulk  flux  calculations  with 
neutral  values  of  0.0014  and  0.0011  for  the  latent  and 
sensible  heat  fluxes,  respectively.  The  evaporative  moisture 
flux  was  derived  from  the  bulk-calculated  latent  heat  tlux. 
The  use  of  bulk  formulas  and  the  model  SST  to  compute  the 
latent  and  sensible  heat  fluxes  provides  an  indirect  relaxa¬ 
tion  to  the  SST  analysis  used  for  GOA M PS.  In  addition,  the 
ocean  model  SST  was  relaxed  to  the  COAMPS  SST 
analysis  via  a  heat  tlux  computed  as  the  difference  between 
the  analysed  SST  anti  the  ocean  model  SST  multiplied  by 
(a  rate  of)  I  m/d. 

3.5.  River  and  Runoff  Inflows 

[it]  Rivers  are  input  into  the  ocean  model  as  a  volume 
source  with  a  specified  vertical  distribution,  zero  salinity, 
and  a  specified  temperature.  Since  real-time  river  temper¬ 
atures  were  not  available,  the  river  temperatures  were 
specified  from  the  monthly  climatology  for  the  Adriatic  at 
the  location  of  the  river  mouth.  River  and  runoff  inflows  for 
the  Adriatic  were  taken  from  the  monthly  climatological 
database  of  Raieich  [1994].  This  database  includes  dis¬ 
charges  for  about  39  rivers  and  runoff  inflows  along  a 
number  of  sections  of  the  Adriatic  coastline.  For  input  to  the 
model,  the  runoff  inflows  were  distributed  along  the  appro¬ 
priate  part  of  the  coast.  The  total  annual  mean  discharge  for 
the  Adriatic  from  Raicich’s  database  is  about  5700  nrVs. 
Daily  observed  discharge  values  w'ere  used  only  for  the  Po 
River  (R.  SignelL  personal  communication,  2003).  The  Po 
was  input  at  5  different  locations  with  each  location  getting 
a  fixed  fraction  of  the  total  Po  discharge. 

3.6.  Model  Simulation 

[is]  The  ocean  model  wfas  run  from  1  September  2002  to 
the  end  of  March  2003  and  model  fields  were  saved  every 
3  hrs  for  analysis  and  comparison  with  observations.  Addi¬ 
tionally,  hourly  values  of  the  model  fields  were  saved  at  the 
mooring  locations  (Figure  2).  The  main  focus  of  the  results 
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Table  I.  Model  Tidal  SSH  Mean  and  RMS  Errors  Relative  to  II 10 


Stations 

Tide 

Amplitude,  cm 

Phase,  deg 

Mean 

RMS 

Mean 

RMS 

M2 

-0,9 

1.2 

-5.3 

8.5 

K! 

- 12 

1.5 

5.3 

7.0 

$2 

0.3 

0,7 

- 15. 1 

16.5 

PI 

-1.0 

1.3 

32.0 

35,0 

Ol 

~0,I 

0.5 

- 12  R 

IK  0 

N2 

0.6 

0.9 

17.8 

26.4 

K2 

-0,3 

04 

42,5 

42.9 

01 

n.i 

0,1 

“3.9 

8.1 

presented  in  this  paper  is  the  time  period  from  the  beginning 
of  January  to  the  end  of  February  2003. 

4*  Observations 

fi<?]  There  were  a  large  number  of  observations  of  differ¬ 
ent  kinds  taken  during  this  experiment,  and  many  of  these 
are  discussed  in  other  papers  in  this  special  issue.  In  this 
paper,  the  ocean  model  results  are  mainly  compared  with 
some  of  the  ADCP  and  CTD  measurements, 

[20]  A  number  of  ADCP  moorings  were  deployed  under 
the  NRL  Adriatic  Circulation  Experiment  (ACE)  together 
with  the  NATO  Undersea  Research  Centre  (NURC)  as  a 
Joint  Research  Project  (JRP),  The  locations  of  these  moor¬ 
ings  are  shown  in  Figure  2+  There  were  14  JRP  moorings 
and  an  additional  upward  looking  ADCP  was  mounted  near 
the  base  of  a  meteorological  tower  [ Cavaleri ,  2000].  The 
J  R  P  m  00  ri  ngs  con  s  i  sted  o  f  t  raw  1  -  re  s  i  s  tan  t ,  bottom  -  moun  ted 
ADCPs  [Perkins  et  ai ,  2000],  All  of  these  moorings 
measured  bottom  temperature  and  pressure  and  three  of 
the  moorings  also  measured  bottom  salinity.  The  ADCPs 
measured  currents  within  the  water  column  from  a  distance 


Figure  3.  Model  and  ADCP  M2  tidal  ellipses  at  10  m. 


Figure  4.  Model  and  ADCP  K1  tidal  ellipses  at  10  m. 


of  about  3  m  above  the  bottom  to  about  5  m  below  the 
surface. 

[21]  The  currents  observed  at  VR5  were  in  fairly  good 
agreement  with  the  model  currents  except  for  the  direction. 
The  disagreement  in  the  direction  was  especially  noticeable 
in  the  tidal  ellipses,  which  were  rotated  about  30°  counter¬ 
clockwise  relative  to  the  tidal  ellipses  predicted  by  NCOM 
and  other  models  (M.  Kuzmic,  personal  communication, 
2004)  and  derived  from  drifters  (P.  Poulain,  personal 
communication,  2005).  Investigation  of  the  VR5  mooring 
data  did  not  reveal  a  reason  for  the  discrepancy.  Because  of 
these  direction  discrepancies,  the  VR5  currents  were  rotated 
30°  clockwise  for  all  the  comparisons  performed  in  this 
paper. 

[22]  The  CTD  data  were  collected  during  several  cruises 
funded  by  the  ADR  1  COS  M,  Euro  ST  RATA  FORM,  and 


Table  2.  Percent  of  Total  Current  Variance  Due  to  Tides  at 
Moonng  Locations 

Name 

5-rrt  Depth 

Near  Bottom 

Model 

ADCP 

Model 

ADCP 

SS2 

17  2 

16.8 

25.5 

20.6 

SS4 

16.8 

15.1 

319 

18.7 

SS5 

21,5 

18.3 

35.7 

29.4 

SS6 

2H.I 

23.0 

43.5 

36.4 

SSH 

30.5 

24,5 

46,7 

39.0 

SS9 

33.2 

24.8 

46  0 

32.4 

ssm 

36.8 

31.2 

40.4 

34.6 

KBI 

5.H 

9.7 

15. 1 

18,2 

CP2 

17.6 

17.5 

40.5 

34,6 

CP3 

25.3 

29.5 

47.4 

42.0 

VRI 

9.7 

10,1 

11.6 

13  7 

VR2 

12.3 

9.7 

20,0 

13,3 

VR4 

22.7 

20.1 

37.0 

35,7 

VR5 

38,5 

30.1 

519 

39.4 

VR6 

61.5 

53,4 

69.7 

61.2 
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Table  3.  Model  Mean,  RMS.  and  Correlation  Errors  With  Respect 
to  ADCP  Measured  Velocities  at  10  ma 


Principal  Axis _  _ Minor  Axis 


Name 

Mean 

RMS 

Cor 

Mean 

RMS 

Cor 

SS2 

4.9 

13,0 

0.601 

-G.t 

5.5 

-0.088 

SS4 

0.1 

8.2 

0,81 1 

-0,1 

5,4 

0  031 

SS5 

2.0 

8.3 

0.732 

-0.9 

5.2 

0.065 

SS6 

5,0 

8.8 

0.683 

-3.9 

7.8 

-0.114 

SS8 

0,4 

5.9 

0.694 

-3.5 

7.4 

0.217 

S39 

-1.4 

5,9 

0,665 

-4.5 

7,6 

0.329 

SSIO 

4,1 

8,4 

0,533 

—0.5 

4,7 

0448 

KBI 

-16.2 

23.2 

0.248 

3,0 

9,6 

OJII 

CP2 

-3.4 

9.3 

0.164 

L5 

6.7 

0.051 

CP3 

-3.5 

9.0 

0.499 

-1.8 

5.8 

0.153 

VRI 

-4.0 

9.2 

0.654 

-0.3 

4.3 

0.034 

VR2 

-0.6 

7.3 

0.645 

-0.2 

3.6 

0.100 

VR4 

1.9 

10.3 

0288 

-2.5 

6.0 

0.094 

VR5 

-0.9 

7.3 

0.621 

-0,8 

5,7 

0.059 

VR6 

0,6 

7.6 

0.275 

2.3 

4,0 

0.475 

“Velocities  arc  dettded  Errors  lire  computed  along  the  local  principal  and 
minor  variance  axes.  Mean  and  RMS  errors  are  in  em/s. 


ACE  projects.  The  locations  of  the  CTDs  taken  during 
January -February  2003  in  the  northern  Adriatic  are  shown 
in  Figure  2, 

5.  Results 

5 A.  Tidal  Sea  Surface  Height 

[23]  The  model  tidal  SSH  was  compared  with  data  from 
27  International  Hydrographic  Organization  (11 IO)  stations. 
These  included  all  the  IHO  stations  in  the  Adriatic  except 
those  in  Venice  Lagoon  (Venice  Lagoon  is  not  adequately 
resolved  by  the  model  grid  and  the  tides  are  not  properly 
simulated  there).  A  harmonic  analysis  of  the  model  tides 


Figure  5,  Derided  model  and  ADCP  mean  current  vectors 
and  std  current  ellipses  at  5  m. 


Figure  6.  Derided  model  and  ADCP  mean  current  vectors 
and  std  current  ellipses  at  30  m. 


was  performed  using  least  squares  titling  to  the  K  tidal 
constituents  used  for  the  tidal  forcing, 

[24]  The  tidal  SSH  mean  and  root  mean  square  (rtns) 
errors  with  respect  to  the  IHO  stations  are  listed  in  Table  1 . 


Figure  7*  Model  mean  currents  for  January  February  at 
5  m. 
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Figure  8*  Detided  model  and  ADCP  10-m  currents  at  SS4  for  January  February  along  axes  of 
(a)  maximum  and  (b)  minimum  variance. 


For  the  largest  tide,  M2,  the  mean  and  rms  errors  for  the 
amplitude  were  —0.9  and  1.2  cm,  respect ively,  and  the 
mean  and  RMS  errors  for  the  phase  were  -53°  and  8,5°. 
The  large  size  of  the  mean  errors  relative  to  the  RMS  errors 
tor  some  of  the  constituents  suggests  that  a  uniform  adjust¬ 
ment  of  the  amplitude  and  phase  of  the  tidal  boundary 
forcing  could  significantly  reduce  the  differences. 

[25]  Several  single-layer,  barotropic  tidal  simulations 
were  performed  on  the  I -km  grid  to  investigate  the  relative 
importance  of  the  tidal  potential  forcing.  Turning  off  the 
tidal  potential  reduced  the  M2  tidal  amplitude  in  the 
northern  Adriatic  about  12%,  e.g.,  the  M2  amplitude  at 
Trieste  decreased  from  26.1  to  23.0  cm.  Turning  off  the  tidal 
potential  reduced  the  K I  and  S2  tides  by  about  7  and  8%, 
respectively. 

5.2.  Tidal  Currents 

[26]  Tidal  currents  for  the  main  tidal  constituents  were 
computed  from  the  mooring  data  and  from  the  model  results 
at  the  mooring  locations  using  least  squares  fitting.  Figures  3 
and  4  show  plots  of  the  model  and  mooring  M2  and  K1  tidal 
current  ellipses.  There  is  fairly  good  agreement  in  the 
amplitude  and  orientation  of  the  tidal  ellipses  at  all 
the  mooring  locations  that  were  used  (the  direction  of  the 
currents  at  VR5  was  modified  as  discussed  in  Section  4). 


The  tidal  currents  for  the  Kl  and  01  constituents  also  agree 
well  with  the  mooring  currents.  The  fraction  of  the  total 
variance  of  the  ADCP  currents  due  to  the  tides  (Table  2) 
ranges  from  a  lowf  of  10%  near  the  surface  at  KB  I  to  64% 
near  the  bottom  at  VR6. 

5.3.  Detided  Currents 

[27]  Table  3  lists  error  statistics  for  the  difference  between 
the  model -predicted  and  ADCP  currents  at  the  moorings  for 
January -February'  at  10-m  depth.  For  all  the  results  pre¬ 
sented  in  this  section,  the  model  and  ADCP  currents  were 
detided  by  subtracting  the  least  squares  fit  tidal  currents. 

|2r]  Figures  5  and  6  show  plots  of  the  5-m  and  30-m 
mean  current  vectors  and  standard  deviation  (std)  ellipses 
for  the  model  and  ADCP  currents  at  the  mooring  locations 
for  January -February,  Note  that  the  scaling  of  the  std 
ellipses  is  smaller  than  for  the  plots  of  the  tidal  ellipses  in 
Figures  3  and  4  by  a  factor  of  2.5. 

[iv]  The  mean  currents  in  the  alongshore  WAC  off 
Senigallia  and  Venice  Lagoon  are  in  good  agreement  with 
the  ADCP  data.  There  are  significant  discrepancies  at  some 
of  the  other  locations.  The  largest  discrepancy  is  at  KIM 
where  the  model  shows  a  much  larger  mean  flow  out  of 
Kvamer  Bay  than  was  observed.  The  model  currents  at  KIM 
are  generally  out  of  the  bay,  whereas  the  observed  currents 
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Figure  9,  Detided  model  and  ADCP  10-m  currents  at  KB1  for  January  February  along  axes  of 
(a)  maximum  and  (b)  minimum  variance. 


are  more  evenly  divided  between  outflowing  and  inflowing 
currents, 

[w]  The  model  currents  at  SS8  and  SS9  just  south  of 
Kvamer  Bay  are  more  directed  toward  the  mouth  of  the  Bay 
than  the  observed  currents, 

[jl]  The  model  mean  currents  are  in  good  agreement  with 
the  observed  mean  currents  at  VR5  in  both  magnitude  and 
direction,  which  supports  the  decision  to  correct  the  direc¬ 
tion  of  the  VR5  currents, 

[32]  The  magnitude  and  orientation  of  the  std  ellipses  of 
the  model  currents  are  generally  in  good  agreement  with 
those  of  the  observed  currents.  The  largest  discrepancy  is 
again  at  KB  l  where  the  std  of  the  model  currents  is  about 
50%  larger  than  that  of  the  observed  currents, 

[n]  Figure  7  shows  the  model  mean  currents  for 
January  - February  at  10-m  depth.  The  mean  circulation 
during  this  period  mainly  reflects  the  model's  mean 
response  to  bora  with  two  large  cyclonic  cells  in  the 
northern  Adriatic  separated  by  a  wedge  of  weaker  cur¬ 
rents  near  1  stria. 

[34]  Figure  8  shows  a  comparison  of  the  10-m  model 
and  ADCP  velocities  at  SS4  for  January- February,  The 
currents  are  plotted  as  components  along  the  axes  of 
maximum  and  minimum  variance,  which  for  SS4  (and 
many  of  the  other  mooring  locations)  correspond  roughly 


to  the  alongshore  and  cross -shore  directions,  respectively, 
as  shown  in  Figures  5  and  6.  The  model  captures  the 
longer-timeseale  (multiday)  events  in  the  WAC  fairly 
well.  These  events  are  bora- induced  intensifications  of 
the  WAC  similar  to  those  described  by  Book  ei  a  I 
[2005]  for  the  winter  of  2001.  However,  there  is  less 
correspondence  between  short-timescale  events  in  the 
model  and  the  ADCP  observations  of  the  WAC.  Along 
the  axis  of  minimum  variance,  the  magnitude  and  time- 
scale  of  the  current  fluctuations  tend  to  be  small  and 
Table  3  indicates  that  there  is  little  correlation  between 
the  model  and  ADCP  currents  along  this  axis  at  many  of 
the  moorings, 

[35]  Figure  9  shows  a  comparison  of  the  10-m  model 
and  ADCP  velocities  at  KB  I.  The  model  shows  outflow 
from  Kvamer  Bay  most  of  the  time  and  strong  outflow 
during  bora  events.  In  contrast,  the  observed  current  is 
often  directed  into  Kvamer  Bay  and  tends  to  be  weaker 
during  boras.  The  model  response  for  simulations  that 
were  run  on  different  grids  with  slightly  different  ba- 
thy  metrics  (section  3.1),  vvas  similar  to  the  model  re¬ 
sponse  shown  here,  which  suggests  that  the  problem  is 
not  due  to  the  specific  geometry  of  the  ocean  model  grid 
that  was  used.  It  may  that  that  there  is  a  problem  with  the 
location  and/or  structure  of  the  bora  wind  jets  over 
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Figure  10,  Model  and  CTD  mean  (a)  temperature  and  (b)  salinity  profiles  for  January  February. 


Kvamer  Bay  predicted  by  CO  AM  PS.  This  discrepancy 
needs  further  investigation. 

5.4,  Comparison  With  CTD  Observations 

[36]  A  large  number  of  CTD  profiles  were  taken  in  the 
northern  Adriatic  during  ACE.  Figure  10  shows  the  mean 
and  Figure  1 1  the  std  of  345  T  and  S  profiles  from  the  CTD 
observations  (Figure  2)  and  from  the  model  simulation 
during  January- February.  For  this  comparison,  the  model 
T  and  S  fields  were  interpolated  to  the  location  and  time  of 
the  CTDs.  Figure  12  shows  the  RMS  differences  between 
the  model  and  CTD  T  and  S  profiles. 

[37]  The  increase  of  the  mean  temperature  with  depth  for 
both  the  model  and  the  CTD  observations  (Figure  10) 
reflects  the  colder  water  in  the  shallower  regions  as  well 
as  temperature  inversions  caused  by  surface  cooling  and  the 
stabilizing  effect  of  the  low-salinity  water  from  the  rivers 
near  the  surface  in  some  areas.  There  is  a  cold  bias  of  0,5  to 
l°C  in  the  model  temperatures  below  20  m.  The  model 
captures  the  mean  increase  of  salinity  with  depth  in  the 
upper  15  m  quite  well,  but  the  model  has  a  low-salinity  bias 
of  0.2  to  0,4  psu  at  depths  below  5  m. 

[3s]  The  std  of  the  model  temperature  profiles  is  similar 
to  that  of  the  CTDs  below  30  m,  but  between  30  m  and 
the  surface  the  difference  steadily  increases,  with  the 
model  std  being  lower.  The  lower  std  of  the  model 
temperature  near  the  surface  is  probably  at  least  partly 
due  to  the  relaxation  of  the  model  SST  to  the  COAMPS 
SSI  analysis.  The  SST  analysis  was  fairly  smooth  and 


hence  tended  to  suppress  variability  in  the  model  near- 
surface  temperature. 

[34]  The  std  of  the  model  salinity  profiles  agrees  very 
well  with  that  of  the  CTDs  (Figure  11).  This  and  the  good 
agreement  of  the  shape  of  the  mean  salinity  profiles 
(Figure  10)  suggests  that  the  model  freshwater  inflows  into 
the  Adriatic  result  in  a  realistic  vertical  salinity  structure, 
with  most  of  the  variability  of  salinity  due  to  the  freshwater 
discharges  being  contained  within  the  upper  15  m. 

[40]  The  small  RMS  error  of  the  model  salinity 
(Figure  12)  relative  to  the  salinity  std  (Figure  11)  and  the 
high  index  of  agreement  [  Willmott  el  ai ,  1985]  of  0,85  for  the 
model  and  observed  salinity  in  the  upper  10  tn  indicate  that 
the  model  simulation  shows  skill  in  predicting  the  horizonal 
variability  of  salinity  structure  due  to  the  freshwater  inflows. 

[41]  The  model  RMS  error  relative  to  the  CTDs  is  about 
l°C  for  temperature  and  about  0.4  psu  for  salinity  below 
15  m.  A  large  part  of  this  RMS  error  is  due  to  the  model's 
low-temperature  and  low-salinity  biases, 

5,5.  Model  Heat  Budget  for  Northern  Adriatic 

[42]  A  heat  budget  for  January  February  was  computed 
for  the  model  tor  the  region  of  the  Adriatic  north  of  the  SS 
mooring  line  (Figure  2).  The  mean  temperature  of  this 
region  decreased  about  4.1°C  during  the  two  months.  The 
effective  temperature  decrease  due  to  surface  heat  fluxes 
was  6.3°C  (which  is  equivalent  to  a  mean  surface  heat  llux 
of  about  175  W/nr),  and  this  was  partially  offset  by  a  mean 
increase  in  temperature  due  to  advection  and  river  inflows 
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Figure  II*  Model  and  CTD  (a)  temperature  and  (b)  salinity  std  profiles  for  January  -  February* 


of  2°C  and  0.2°C,  respectively.  The  use  of  proper  river 
in  How  temperatures  would  likely  have  led  to  a  net  cooling 
from  the  rivers,  rather  than  a  slight  warming. 

5,6,  Bora  Response 

[4.1]  The  effect  of  the  winds  and  wind  stress  curl  during 
bora  events  has  been  noted  by  a  number  of  investigators 
[Or lie  ei  at,  1994;  Rachev  and  Purini ,  2001 ;  Beg  Paklar  et 
ai,  2001].  The  variations  in  the  mountain  topography  of 
Croatia  act  to  force  the  bora  winds  out  over  the  water  in 
strong  jets  and  set  up  a  wind  stress  pattern  that  tends  to 
recur  during  different  bora  events  with  variations  mainly  in 
amplitude  and  duration. 

1 44 1  The  COAMPS  winds  used  for  the  model  simulations 
here  have  sufficient  resolution  (4  km)  to  respond  to  the  main 
features  of  the  orography  along  the  Croatian  coast  and 
provide  a  fairly  realistic  pattern  of  strong,  separated  jets 
[Pullen  et  ai,  2003],  Figure  13  shows  a  plot  of  the 
COAMPS  24-hour-averaged  wind  stress  during  the  bora 
of  12  January. 

[4*]  Figure  14  shows  the  model  response  to  the  bora  that 
occurred  on  12  January.  The  response  to  this  bora  by  the 
ocean  model  simulations  run  on  other  computational  grids 
was  very  similar,  and  other  bora  events  during  the  winter  of 
2003  generated  similar  circulation  patterns. 

[46]  Figure  14  shows  large,  cyclonic  gyres  in  the  north- 
central  and  nonhem  Adriatic  and  a  smaller  ant  [cyclonic 


gyre  near  the  Istrian  Peninsula  between  the  large  cyclonic 
gyres.  This  double-gyre  pattern  is  consistent  with  past 
theoretical  and  modeling  studies  and  the  JRP  data  offer 
clear  evidence  for  the  existence  of  the  two  strong  cyclonic 
cells. 

[47]  The  model  response  to  the  bora  of  12  January 
produces  circulations  in  the  southern  half  of  Figure  14  that 
depart  significantly  from  the  classic  view  of  a  large-scale 
cyclonic  circulation  trained  by  the  WAC  and  the  EAC.  The 
WAC  is  present  ail  along  the  Italian  coast,  hut  a  significant 
portion  of  its  southeastward  tlow  north  of  Ancona  recircu¬ 
lates  toward  Croatia  forming  a  tighter  cyclone.  Below'  this 
cyclone,  the  EAC  has  been  replaced  by  a  small  cyclone  - 
anticyclone  pair.  Hence,  during  this  time  period  of  intense 
wind,  there  are  five  separate  semiciosed  cells  in  the  northern 
Adriatic. 

[4*]  The  closing  of  the  two  most  southern  cells  is 
transient.  The  less  frequent  of  these  two  cells  is  the  more 
southerly,  cyclonic  cell.  However,  this  cell  appears  in  the 
model  simulation  at  least  once  (though  usually  briefly) 
during  each  of  the  five  bora  that  occurred  during  January  - 
February.  When  these  gyres  are  not  present  during  bora 
winds,  there  are  frequently  corresponding  cyclonic  and 
anticyc Ionic  meanders  in  the  EAC  and  these  meanders 
are  persistent  enough  to  be  preserved  in  the  January' 
February  mean  field  (Figure  7), 
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Figure  12.  Model  and  CTD  RMS  (a)  temperature  and  fb)  salinity  profile  differences  for  January 
February. 


[49]  Unfortunately,  the  measurements  taken  during  this 
period  were  not  located  in  optimal  positions  for  validating 
the  tight  recirculation,  the  meanders  of  the  EAC,  and  the 
closing  of  these  meanders  to  form  the  separate  cyclonic  and 
anticyclonic  cells  that  are  observed  in  the  model  simula¬ 
tions.  However,  the  data  do  present  some  evidence  in 
support  of  these  features.  AL  the  start  of  each  strong  bora 
period,  the  observed  barotropic  velocities  at  mooring  SS8 
have  significant  flow  toward  Croatia.  That  is,  velocities  are 
briefly  directed  approximately  toward  5°T  rather  than  to¬ 
ward  — 40°T.  During  bora,  the  shear  in  the  alongshore 
velocity  component  between  moorings  SS10  and  SS9  peaks 
in  an  anticlockwise  sense.  This  result  supports  the  model- 
simulated  anticyclonic  cell  or  meander  in  the  EAC  south  of 
these  positions.  ADCP  measurements  made  by  the  EACE 
program  [Orlic  ei  ai ,  2006]  show  cyclonic-type  shear  in  the 
EAC  for  January  and  February.  The  southernmost  cyclonic 
cell  in  Figure  14  is  positioned  at  the  sites  of  these  two 
moorings,  and  thus  this  shear  lends  support  to  the  predic¬ 
tions  of  a  cyclonic  meander  or  cell.  However,  peaks  in  shear 
between  these  positions  from  the  data  or  from  the  model  do 
not  showr  the  same  level  of  correlation  with  bora  events  as 
the  results  from  SS10  and  S$9. 

[so]  Figure  15  shows  a  plot  of  the  COAMPS  wind  stress 
curl  averaged  over  the  24  h  proceeding  the  circulation 
snapshot  shown  in  Figure  14.  The  areas  of  positive  and 
negative  wind  stress  curl  match  the  regions  of  cyclonic  and 


anticyclonic  circulation  in  Figure  15  and  show-  the  strong 
control  exerted  by  the  winds  on  the  circulation  during  these 
events, 

6.  Summary  and  Conclusions 

[51]  The  NCOM-modeled  tidal  height  shows  good  agree¬ 
ment  with  IHO  stations  in  the  Adriatic  and  the  modeled  tidal 
current  ellipses  show-  good  agreement  w  ith  ADCP  measure¬ 
ments,  The  fraction  of  the  total  current  variance  due  to  the 
tides  at  the  mooring  locations  at  5-m  depth  ranges  from  6  to 
70%. 

[5:]  Comparison  of  derided  model  and  ADCP  currents  at 
10-m  depth  for  January  -  February  show's  RMS  errors  along 
the  principal  variance  axes  ranging  from  6  to  12.9  cm/s  and 
correlations  ranging  from  0.16  to  0.81.  Correlations  be¬ 
tween  the  detided  mode!  and  ADCP  currents  along  the 
minor  variance  axes  are  generally  low. 

[53]  In  general,  the  model  predicts  longer-timescale  fea¬ 
tures  well  at  all  the  JRP  locations  except  at  KB1.  The  lack 
of  model-to-data  correlation  between  shorter-timescale  fea¬ 
tures  and  minor  variance  axis  flows  may  be  due  to  the 
presence  of  small-scale  eddies  and  instabilities.  The  phasing 
and  location  of  such  features  may  not  be  very  deterministic 
as  demonstrated  by  the  sensitivity  of  these  features  to  the 
different  ocean  model  grids  used  in  this  study.  The  dis¬ 
agreement  between  the  model  and  the  measurements  at  KB  I 
requires  further  study  and  explanation. 
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Figure  13,  Time-averaged  COAMPS  wind  stress  for 
period  09Z,  11  January,  to  09Z,  12  January  2003.  Contour 
interval  is  0J  Pa, 


Figure  14,  Model  IG-m  currents  at  09Z,  12  January  2003, 


Figure  15.  Time-averaged  CO  AMPS  wind  stress  curl  tor 
period  09Z,  1 1  January,  to  G9Z,  12  January  2003.  Contour 
interval  is  0,01  Pa/km. 


[54]  Comparison  of  the  model-simulated  temperature  and 
salinity  profiles  during  January -February  with  CTD  mea¬ 
surements  indicate  that  the  model  captures  some  of  the 
spatial  structure  of  the  observed  fields.  The  good  agreement 
of  the  model  and  observed  salinity  mean  and  std  and  the 
high  (0.85)  index  of  agreement  of  the  model  salinity  with 
the  observed  salinity  in  the  upper  10  m  indicate  that  the 
model  freshwater  inflows  into  the  Adriatic  result  in  a 
realistic  vertical  salinity  structure  and  that  the  model  has 
some  skill  in  predicting  horizontal  salinity  variations.  The 
model  show  s  a  bias  of  about  —  1  aC  below'  20  m  and  -0,2  to 
—0.4  psu  relative  to  the  CTDs. 

[55]  The  model  response  to  the  several  bora  events  in 
January'  February  shows  a  recurring  pattern  of  large  cyclo¬ 
nic  and  smaller  anti  eye  Ionic  gyres  that  generally  agree  with 
the  ADCP  observations  and  reflect  the  pattern  of  w  ind  stress 
curl  from  the  COAMPS  wind  stress  forcing.  During  strong 
bora  forcing,  two  large  cyclonic  circulation  gyres  form  in 
the  northern  Adriatic  with  a  smaller  anticyclonic  circulation 
between  them  near  the  Istrian  Peninsula,  Southeast  of 
K varner  Bay,  a  couple  of  relatively  large  meanders  fre¬ 
quently  occur  within  the  Eastern  Adriatic  Current  (EAC) 
and  these  meanders  sometimes  dose  to  form  small  cells. 
The  circulation  within  the  meander/cell  closer  to  K varner 
Bay  is  anticyclonic  and  the  more  southeasterly  meander/cell 
is  cyclonic. 
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